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ABSTRACT 

The shocked interstellar medium around IC443 produces strong X-ray emission in the soft energy 
band (E < 1.5 keV). We present an analysis of such emission as observed with the EPIC MOS cameras 
on board the XMM-Newton observatory, with the purpose to find clear signatures of the interactions 
with the interstellar medium (ISM) in the X-ray band, which may complement results obtained in 
other wavelenghts. 

We found that the giant molecular cloud mapped in CO emission is located in the foreground and 
gives an evident signature in the absorption of X-rays. This cloud may have a torus shape and the 
part of torus interacting with the IC443 shock gives rise to 2MASS-K emission in the southeast. The 
measured density of emitting X-ray shocked plasma increases toward the northeastern limb, where 
the remnant is interacting with an atomic cloud. We found an excellent correlation between emission 
in the 0.3-0.5 keV band and bright optical/radio filament on large spatial scales. The partial shell 
structure seen in this band therefore traces the encounter with the atomic cloud. 
Subject headings: X-rays: ISM — ISM: supernova remnants — ISM: individual object: IC 443 



1. INTRODUCTION 

SNR IC443 (G189.1-h3.0) is located in the Gem OBI 
assoc iation at a distance of 1.5 kpc ijWelsh fc SallmenI 
UnnS)- Optical and radio morphology is shell- like, it 
has an angular diameter of 50' and appears to con- 
sist of two connecte d subshells with different radii 
()Braun fc Stromlll986() . In the X-ray band IC 443 has 
been previously observed with Ein stein and the HEAO 1 
A2 experiment jPetre et allllOSSD, Ginga llWang etaP 
[1992'), ROSAT fAsaoka & Aschenbach" '1994J), ASCA 
iKcolianc ct al. 19 97; Kawasaki ct al. 2002) , BeppoSAX 
jBocchino fc Bvkovl 1200(1) and RXTE (|Sturner et all 
l200ir 

IC443 shows features simila r to the class of mixe d mor- 
phology supernova remnant ijRho fc Petrel 1199^ . The 
X-ray morphology is centrally peaked and there is no ev- 
idence of a limb-brightened X-ray shell. The X-ray emis- 
sion is primarily thermal and it has been described with 
a two component ionization equilibrium model of tem- 
peratures ~0.2-0.3 keV (hereafter cold component) and 
~1.0 keV (hereafter hot co mponent, Pctre ct al. 1988; 
lAsaoka fc Asch enbach"1994l. On the contrary, the work 
of lKawasakTet al. (2002) suggests a non-equilibrium con- 
dition of the X-ray emitting plasma, reporting a low ion- 
ization timescale for the cold component and a ovcrion- 
ization state for the hot one. 

Recently Chandra and XMM-Newton observations 
identified a plerion nebula in the southern part of th e 
remnant (lOlbert et al.ll200Tl IBocchino fc Bvkovl l200T]). 
but its association with IC443 is still under debate. 
Source position seems not compatible at the 90% con- 
fidence limit w ith the EGRET sourc e 3EG J0617-h2238 
llEsDOsito et all. 1996.) ■ as argued bv IBocchino fc Bvkovl 

^ Dipartimento di Scienze Fisiche ed Astronomiche, Sezionc di 
Astronomia, Universita di Palermo, Piazza del Parlamento 1, 90134 
Palermo, Italy 

^ INAF-Osservatorio Astronomico di Palermo, Piazza del Parla- 
mento 1, 90134 Palermo, Italy 



The environment in which the remnant is evolv- 
ing seems to be quite complex. In the south the 
blast wave has been decelerated by the encounter with 
a dense (nn^lO"^ c m~^) and clumpy m olecular ring 
(jBurton et a l.lll988l: IPickman et alJll992t). A rich spec- 
truni of mo lecula r lin es emission l)van Dishoeck et alJ 
119931 iRho et al.l 120011 and references therein) together 
with OH(1720 MHz) mas er emission (Claussen et al] 
Il997t iHoffman et al.l I2003D have been detected, con- 
firming the interaction with high density gas. In the 
northeast the shock front has encountered a less dense 
(10<no<103 cm-3) Hi cloud l|Denoveilll97l IRho et all 

m^- 

Emission from quiescent molecular gas has been observed 
toward IC443 direction and it is likely due to a giant 
molecular cloud in front of the remnant l)Cornett et alJ 
Il977j) . A superposition between IC443 and another 
separated S NR (G189.6+3.3) has also been suggested 
ijAsaoka fc Aschenb ach 1994). 

Multiwavelenght observations have shown the presence of 
sharp density gradients and different clouds geometries 
in the surroundings of IC443. Establishing the exact po- 
sition of ISM structures respect to the IC443 is of crucial 
importance to understand SNR dynamics and evolution. 

In this paper, we focus on the interaction between the 
SNR shock and the dense ISM environment. For this 
reason, we present a detailed analysis of the softer X- 
ray thermal emission of the SNR IC443 based on XMM- 
Newton observatory data. We leave the study of the hot 
thermal component, which is more prominent in the in- 
ner regions of the remnant, to a successive paper, now 
in preparation. The aim of this work is to take advan- 
tage of XMM-Newton characteristics to derive a better 
description of the SNR evolutionary scenario and to di- 
rectly measure physical conditions of the X-ray emitting 
plasma both in the bright regions and in the yet unex- 
plored faint western regions. 

In Section 13 we present our observations and the data 
reduction methods we used. Section 12.21 provides a de- 
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TABLE 1 

General informations about MOS data (MOS1+MOS2). 



OBSJD 


Date 


"J2000 


'5j2000 


^ exp 






UT 


(hh:mm:ss) 


(dd:mm:ss) 


(ksec) 




0114100101 


26/09/2000 


06:17:28 


+22:41:44 


20/45 




0114100201 


25/09/2000 


06:16:15 


+22:41:60 


10/11 


■a- 


0114100501 


25/09/2000 


06:16:15 


+22:41:60 


40/50 


1 


0114100601 


27/09/2000 


06:17:28 


+22:25:44 


11/12 




0114100301 


27/09/2000 


06:17:28 


+22:25:44 


40/50 




0114100401 


28/09/2000 


06:16:15 


+22:18:00 


50/60 


i 



^Screened/unscreened exposure time 

scription of our background subtraction method. In Sec- 
tion O we present the results obtained from the XMM- 
Newton data analysis. Section^contains a discussion on 
the physical properties of the cold X-ray emitting plasma 
fSect. H!T|l . on the interaction between the remnant and 
the surrounding clouds fSect. lO|) and on the absorption 
from a foreground dense molecular cloud (Sect. In 
Section [SJ we summarize our results and conclusions. 

2. DATA ANALYSIS 
2.1. Observations 

IC443 was observed on 2000 September as part of the 
Cal/PV phase of the XMM- Newton Observatory, using 
EPIC and RGS instruments. The observation was per- 
formed in four distinct pointings to allow a maximum 
coverage of the source. These are the same observations 
analyzed bv lBocchino fc Bvkovl l|2003|) . 

Here, we focu s only on the data c ollected by the EPIC- 
MOS cameras l|Turner et alJl200]]) . During each obser- 
vation, the MOS cameras were operated in Full Frame 
Mode, providing a temporal resolution of 2.6 s, and the 
medium filter was used to limit the number of low-energy 
photons. The Science Analysis System software (SAS, 
version 6.0) was used for data reductio n. Spectral an al- 
ysis was performed with XSPEC v. 11.1 ijArnau 

We generated calibrated events files with the SAS task 
emproc and then we further screened the data, selecting 
only events with PATTERN < 12 and an optimal value of 
FLAG, as suggested by the current status of calibration'^. 
We discarded events with energy below 0.3 keV from our 
future analysis because of data calibration uncertainties. 

2.2. Background subtraction 

The EPIC MOS background components are mainly 
electronic noise, fluorescence emission lines (instrumen- 
tal background, not vignetted). X-ray galactic emission 
(astrophysical background, vignetted) and soft magneto- 
spheric protons fo cused by the optics ijLumb et alJl2002t 
IMartv et al.1l2003|) 

We removed high proton background time intervals by 
extracting light curves at energies > 8 keV, where the 
EPIC cameras have a low quantum efficiency, with a bin 
time of 30 s . We rejected data intervals with more than 
15 counts. 

To check for a residual proton component we 
used a simple but usefu l diagnostic, developed by 
iDe Luca fc Molendil l)2004j) . We calculated for each 




I 
I 




Fig. 1. — X-ray count rate image of IC443 in the 0.3-1.4 keV en- 
ergy band. The image was adaptively smoothed using crmi„=5"and 
"■mai=20"; the bin size is 5". Contour levels correspond to 1.4, 4, 
7, 10 10"* ctss-i bin-i. Arrows indicate regions of low brightness 
surface; the partial shell structures are labeled as filNE and filNW. 

screened data set the ratio R between the surface bright- 
ness in FOV (S/jv) and the surface brightness out FOV 
CSqut) in the energy range 8-12 keV. Fol lowing the defi- 
nition given in lDe Luca fc Molendil l)2004j) , the maximum 
value of R, above which there is a significant soft proton 
contamination, is R„ia2; = l-3. If 1.05<R<1.3 the residual 
proton level is low, if R<1.05 it can be considered neg- 
ligible. We verified that our cleaned data sets have all a 
low or negligible residual proton component. After the 
screening process, the total exposure time was reduced 
by 25% (from 230 to 170 ks). 

Table ^ summarizes the main information about the 
observation: the date, the pointing's center coordinates, 
the exposure time before and after flare screening. 

Instrumental background is the dominant component 
of the events out of FOV and it is also highly spatial 
variable. A typical solution is to estimate it using merged 
data sets acquired with the filter wheel in the CLOSED 
position ((Lumb et al. 2002) . We used the closed data 
sets of P. Marty, distributed by the School of Physics 
and Astronomy of the University of Birmingham^. 

Since IC443 is an extended source, that fills the en- 
tire field of view, it is not possible to estimate the local 
astrophysical background from the same data. Instead 
of a local background, we used a set of observations of 
blank sky regions, produced from the Galactic Plane Sur- 
vey observations (GPS, P.I. Parmar), pointed toward a 
galactic longitude of 319° and a galactic latitude between 
0° and 3°. This data set represents a better estimate of 
the cosmic background at low galactic latitude than the 
standard blank sky field at high 

After the same screening process described in Sect.|5]and 
the removal of point sources, data sets were scaled by the 
relative exposure time and then reprojected onto the sky 
attitude of IC443 pointings with the script skycast. 



^ http:/ /xmm.esac.csa.int/external/xmm_sw_cal/calib/ 



* http:/ /www. sr. bham.ac.uk/xmm3/scripts. html 
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Fig. 2. — Left panel: X-ray count rate image of IC443 in the 0.3-0.5 keV energy band. The image was adaptively smoothed using 
O'min=^0" and crmai=40"; the bin size is 10". We overlaid VLA NVSS contours at 1.4 GHz, corresponding to 5, 20, 45, 80 mjy beam~^. 
Right panel: Optical DSS image of the northeastern shell of the remnant overlaid with contours of X-ray emission in the 0.3—0.5 keV energy 
band. Contour levels correspond to 1.7 (blue), 2.6 (green), 4.6 (yellow), 8.6 (red) 10""^ cts s"'^ arcmin"^. 

From a comparison between the spectra extracted from 
the whole FOV of the CLOSED, GPS and IC443 data, 
we noted that in the energy range 0.3-5.0 keV the source 
signal is the brightest, while above 5.0 keV the instru- 
mental noise component is the dominant one. 
We decided to consider in our future analysis only events 
in the band 0.3-5.0 keV and to use GPS data as back- 
ground correction. We verified that in this band our 
results are not sensitive to the background subtraction 
(either GPS or high galactic latitude blank fields). 



3. RESULTS 

3.1. X-ray Images 

We defined two energy bands, in which images has been 
extracted, so that in each band the contribution of a sin- 
gle component is the dominant one and, at the same time, 
the statistic data quality is good. In the chosen range 
0.3-1.4 keV (soft band), the cold component fiux (cal- 
culated assumi ng the parameters of the two co mponents 
fit re ported by iPetre et aP 1 19881 lAsaoka fc As chenbacO 
119941) is the 64% of the total; in the range 1.4-5.0 keV 
(hard band), the hot component fiux represents the 90% 
of the total. 

A third image, generated in the 0.3-0.5 keV band, traces 
the spatial distribution of the coldest X-ray emitting 
plasma and it is useful to trace the boundary between 
the SNR shock and the surrounding molecular clouds. 

All the X-ray emission images presented here are back- 
ground subtracted, exposure and vignetting corrected. 
To correct for the exposure and the vignetting, we di- 
vided the mosaiced background-subtracted count image 
by its mosaiced exposure map. The resultant image was 
adaptively smoothed with the SAS task asmooth. 

Figure ^ shows the 4 pointings mosaiced image of 
IC443 in the 0.3-1.4 keV band. The brightest features 



are concentrated in the northeastern (NE) quadrant and 
a clear shell is not present. We mark a region near the 
rim of the remnant (filNE in Figure ^ : the high sur- 
face brightness gradient suggests that we are looking the 
shock front nearly edge-on. The large effective area and 
the good spatial resolution of XMM-Newton allow us to 
have a better view of the western part of the remnant 
and to resolve its structure for the first time. Western 
X-ray emission appears about 10 times fainter and more 
diffuse than the eastern one. We can identify a bright 
filament (filNW in Fig. that suggests the existence of 
a limb-brightened partial shell. 

The arrows in Fig. ^ indicate two dark lanes, that 
cross the remnant: the first, from north to south, is 
probably due to the presence of a molecular clou d be- 
tween the observer and IC443 ijPetre et al.lll988|) : the 
second, fro m northeast to southwest, was explained by 
lAsaoka ~A schenbach ( 1994) as an additional absorption 
effect by the SNR G189.6-H3.3. 

On the other hand, the morphology of the very soft 
emission (0.3-0.5 keV, see Fig. is not centrally filled, 
the X-ray surface brightness peaks are near the north- 
eastern edge of the remnant, resembling strikingly IC443 
radio and optical emission. For comparison, we show in 
Fig.|21(/e/i5 panel) the X-ray image of the remnant in the 
0.3-0.5 keV band, overlaid with contours of radio emis- 
sion lev els at 1.4 G Hz, observed with the VLA NVSS 
survey |Cpiidon_et al. 1998). We also present the su- 
perposition between X-ray contours and t he optical im- 
age, extracted from the DSS red survey l)Lasker et alJ 
119901) . of the northeastern quadrant of IC 443, where it 
is encountering a Hi cloud l)Denoveilll978l) (Fig. El right 
panel). 

In the very soft X-ray band it is possible to resolve a 
partial shell structure in the NE, similar to that we see 
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in the radio and optical images. iKawasaki et all l|2002fl 
had aheady claimed the presence of a shell-like struc- 
ture on the basis of ASCA GIS softness ratio map in 
the energy bands 0.7-1.5 and 1.5-10 keV. However our 
results represent a direct detection of a very soft shell, 
which is likely differen t from the structure reported by 
IKawasaki et al.l 1)2002(1 . In the light of our result on 
molecular cloud obscuration of the IC443 X-ray emis- 
sion (see p ar. 14.31 below), we arg ue that the structure 
detected by IKawasaki et aP f2002") in their softness ratio 
map is due to the absorpt i on of the foreground molecu- 
lar cloud of iCornett et alJ l(1977|) instead of a real shell 
morphology. 

In the northeastern quadrant, radio, optical and very 
soft X-ray emission correlate well: X-ray brightest fea- 
tures are confined within the radio-optical shell and the 
position of maximum X-ray emission contour levels, just 
behind bright optical filaments, is in agreement with a 
shock front that has just encountered a density gradi- 
ent of the ISM with higher density located in the outer 
region s. This correlation was also noted bv iPetre et al.l 
(|1988j) . which however focused only on the northern edge 
of the brightest soft X-ray region. Our results instead 
shows that this correlation involves the very soft X-ray 
shell and the optical-radio emission on a much greater 
extent. 

The structures seen in the 0.3-0.5 keV map are of course 
affected by absorption, however we will see in sect 13. 21 
that the shape of the very soft shell is not correlated 
with the absorption map in the east of the SNR. In the 
light also of the good correlation between the very soft 
shell and the radio emission in the NE of IC443, we argue 
that a real shell structure exists. 

In the southeastern quadrant, where the remnant is in- 
teracting with a dense molecular cloud, the spatial gap 
between X-ray emission and radio contours, as shown 
in Fig. 13 suggests that the shock front has been pen- 
etrated deeply inside the ring of molecular gas and it 
has been stron gly decelerated, em i tting mainly i n the 
infrared band feurton et all Il988t IRichter et all 119951 
IRho et alJl200l . 

3.2. Median Energy Map 

The complicated morphology of IC443 and the previ- 
ous literature results suggest that the plasma has a dis- 
tribution of physical conditions and so it is crucial to 
perform our spectral analysis selecting regions as homo- 
geneous as possible. The median ph oton energy Erq is a 
robust spectral indicator, as shown in lHong et al.l l|2004|) . 
that allows us to map spectral variations inside the field 
of view. 

For a plasma in ionization equilibrium and with solar 
abundances (model MEKAL in XSPEC) we verified that 
in the soft band (0.3-1.4 keV) E50 is mainly influenced 
by absorbing column variations, especially at large values 
of the absorption (Nh>0.5x10^^ cm~^). The choice of 
the band assures that the median energy is not heavily 
affected by the parameters of the X-ray hot component. 

Fig. shows the median energy map of the emitting 
plasma in the 0.3-1.4 keV band. Molecular line data to- 
ward IC443 are plotted as contours over the image: qui- 
escent CO in white (Lee J. J., private communication) 
and shocked II2 in black l|Burton et al.llT988|) . 
To build the map we used a pixel size of 20", in or- 




Fig. 3. — Median energy map in the soft band (0.3—1.4 keV), 
overlaid with quiescent CO (white) and shocked H2 (black) emis- 
sion contours. The arrows indicate two high median energy regions 
without a counterpart in CO or H2 maps. We masked some regions 
inside the observed field of view, traced by the dashed line, as they 
are too noisy (S/N<3). Image bin size is 20" and smoothing width 
is o-=20". The velocity-integrated (-6 to km s"!) J=l-0 CO 
emission has 50" resolution. The velocity-integrated d=1-0 S(1) 
H2 line emission a t 2.122 ^m has the spatial resolution of 40" 
IBurton et aLllTosl) . 




Fig. 4. — Median/rate plot for the soft band (0.3-1.4 keV): on 
the Y-axis is reported the median photon energy, on the X- axis 
the corresponding count rate. The dashed line traces the median 
energy threshold at 0.97 keV between "on cloud" and "off cloud" 
regions. Points come from five regions, superimposed on the me- 
dian energy map in the top right corner. The black contour in the 
map corresponds to the value of E5o=0.97 keV. 

der to collect at least 10 photons per pixel, and then, 
for each pixel we calculated the median photon energy. 
The resulting image was smoothed using a bidimensional 
Gaussian distribution with cr=20" . No correction for the 
X-ray background was made, we simply masked faint re- 
gions, where the signal to noise ratio was less than 3. 
The image presented in Fig.|21allows us to compare spec- 
tral properties of the cold X-ray component with the 
emission of shocked and quiescent molecular gas with an 
unprecedented spatial resolution (28" , corresponding at 

^0.2 pc at a distance of 1.5 kpc). 

Previous studies about IC443 l|Petre et all 119881: 
lAsaoka fc Aschenbachl[T994|) tried to link spectral hard- 
ness ratio variations across the remnant to an absorption 
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effect, but the low images resolution did not permit to 
study the correlation with the molecular cloud material 
on small angular scales. Moreover, previous hardness 
ratio maps are sensitive both to absorption effects and 
to the parameters of the hot component, whereas our 
median energy map is more strongly dependent on ab- 
sorption. 

The spatial variations of the spectral properties reported 
in our map f Fip||5)l confirm the distribution of the hard- 
ness ratio in lAsaoka fc AschenbachI l)1994() . but shows 
more details. In particular, we find a strict correlation, 
also on small angular scales, between molecular emission 
contours and high median energy regions, as expected in 
the case of an increasing column density. They match 
very well everywhere, except two high median regions, 
indicated by the arrows in Fig. 13 that seems not to have 
a counterpart in CO and H2 maps. 

As alrea dy shown in t he stu dy of the Vela supernova 
remnant of IMiceli et al.l l)2005j) . the relation between the 
median photon energy E50 and the corresponding count 
rate R in each pixel is able to give a global view of differ- 
ent thermal structures and spectral variations inside the 
field of view, preliminary to a complete spectral analy- 
sis. As shown in Fig.0] we selected five significant regions 
and made the corresponding E50 vs. R plot to get fur- 
ther information about the physical parameters of the 
plasma. We applied the Maritz-Jarrett method to esti- 
mate the errors o n the median energy, as described by 
iHong et all l|2n04D . and we assumed a poissonian error 
for the count rate values. In order to avoid excessive 
noise, especially in faint areas, points with cr(i?5o) >5% 
or a{R) >15% were not reported in the plot. 

It is worth to note that all the points belonging to re- 
gion 1 (blue circles), which corresponds to a CO emis- 
sion area, have the highest median values (E50 >1.0 
keV), while points from other regions lie below this en- 
ergy threshold. In particular, we can distinguish two 
distinct branches at 1.0 keV (region 1) and at 0.9 keV 
(region 2 and a great part of region 3). These features 
could be interpreted as an absorption effect: on average 
plasma physical conditions appear to be the same in dif- 
ferent regions, the only difference being an increase of the 
column density, that shift points from region 1 towards 
higher median and lower rate values. We may define the 
boundary of off-cloud and on-cloud regions with the value 
E5o=0.97 keV (dashed line in Fig. 0] left panel), whose 
contour is also reported in the right panel of Fig. 01 Data 
from the off-cloud region 4 (yellow diamonds) are those 
affected by the greatest uncertainties. They show similar 
E50 values of region 3 (green squares) , but seem isolated 
from the others, maybe suggesting physical differences, 
e.g. a lower interstellar absorption. The arc-like shape, 
identified by region 5 red points, is characteristic of tem- 

Sierature variations for a plasma in pressure equilibrium 
MiceU et al.l200l . It seems also present in the off-cloud 
region 3, though less evident. 

3.3. Spectral Analysis 

The XMM-Newton data allowed us to significantly im- 
prove the spatial resolution of the spectral analysis with 
respect to past work s about IC443 (Pctrc et al. 1988; 
IKawasaki etalll2002fl . 

We analyzed spectra extracted from the regions marked 
in Fig. [S| These regions were selected so as to represent 
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Fig. 5. — X-ray count rate image of IC443 in the 0.3-1.4 keV 
energy band, on which we superimposed the nine regions selected 
for spectral analysis and the median photon energy contour (white) 
corresponding to E5o=0.97 keV. The dashed semicircles mark sub- 
shell A and subshell B structures. The color bar is the same as in 
Fig-El 

Model: WABS {VPSHOCK + VMEKAL) 




Energy (keV) 

Fig. 6. — XMM-Newton EPIC MOS spectra representative of 
IC443 thermal emission. The top figure shows the spectrum of reg 
1 with the corresponding two components best fit model and, in 
the lower panel, the residuals distribution. The blue dashed line 
represents the cold spectral component, while the red dotted line 
represents the hot one. The bottom figure shows the spectrum of 
reg 8 with its corresponding single component best fit model and, 
in the lower panel, the residuals distribution. 

the main features of the cold plasma, as obtained from 
imaging analysis. We focused on the following aspects: 
(a) to verify the presence of a column density gradient 
across the FOV and its correspondence with median pho- 
ton energy variations (regions 1-5); (b) to characterize 
physically the partial shell structure (regions 7 and 8) 



6 



E. Troja et al. 



TABLE 2 

Dimensions, counts and best-fit parameter values for the nine spectral regions. For each spectrum we also report 

THE derived filling FACTOR FOR EACH SPECTRAL COMPONENT. QUANTITIES IN SQUARE BRACKETS ARE HELD FIXED AT THE GIVEN 

VALUE. All the reported errors are at the 90% CONFIDENCE LEVEL. 

Parameters Reg 1 Reg 2 Reg 3 Reg 4 Reg 5 Reg 6 Reg 7 Reg 8 Reg 9 

Area (arcmin^) . 5.2 1.1 2.2 7.5 2.3 1.1 5.5 5.0 7.1 

Counts (103) ... 5.6 6.9 5.9 5.8 6.3 7.0 4.6 4.8 6.7 

NHdO^^cm-^) 0.96«-i^ 0.74«;}i 0.561"-^ 0.431^1 0.451°;°^ 0.6At°Z 0.57^^11 0.58«;»? 0.6^^ 

Cold component 

kT.(keV) 0.37iHI 0.58«;«^ 0.291^1 O.SltHI OSA^lf, 0.358l°:«f^ 0.291^^ O.SltVol 0.66i«:« 

T(10i2cm-3s) 48 (>2) 2.oti.2 50 {>2) 2.3 (>0.7) CIE 50 {>3.3) 0.2lj;JJ;*o l.oto.7 l-2to.8 

Fs'-'^^ 4^3 6±3 llt^" 8±4 2.3±0.8 7.0^2* 1619" ™-lo 2,1±1,8 

O/O0 [ 1 ] 4.3tl6 l-3±a7 0.44t»:i^ 2.6t? 0.8±» j 0.8t« j 0-52^°% I 1 1 

Ne/Nee 1.9t% 3.5^^ l-2tL 0.451"-? 2.8tJ:^ l-Ot^j O-^tll 0.7^^2 ^-^^U 

Mg/Mg0 l.olg.5 l-Sto.l 0-7lo:4 0.32i0:i? 1.5±0.5 0.48^^2 O-StH 0.50t°;?7 1.41^5 

Si/Si0 [1] [1] [1] [1] 0.4«;?g [1] l.ltJ:^ 0.5l°:| [1] 

S/S© [ 1 ] [ 1 ] [ 1 ] [ 1 ] 0.13 {<0.4) [ 1 ] [ 1 ] [ 1 ] [ 1 ] 

Fe/Fee OA^H 0.9«-j 1.3^^^ 0.32i«:0^ 0.15i°:Jg 0.511°;?, 0.8±0.2 0.451°;?, 0.5^11 

filling factor /s . 0.15 0.42 0.26 1.0 1.0 0.52 1.0 1.0 0.2 

Hot component 

kT^ CkeV~l 1 11+0.3 1 70+0.13 1 r,i+0.i3 ~ ~ 1 4+0.8 ~ ~ 1 70+0.19 

Kl/i (.Key; J^ J^J^-0.16 ^-'^-0.18 ^^■'-'^^-0.06 ^-^-0.3 ^^•'•^-0.13 

p, (a) 1 A+0-8 r, An+0.6 1 c+0.8 p, .r,+0.3 n S+0'5 

Mg/MgQ 2.3l? 1 >5 0.8lg;^ — — >2 — — >4 

Si/Si0 0.5±0.3 >3 1.0±0.5 — — >5 — — l-^lo.e 

S/Sq OJtol >3 l.OitO.6 — — >3 — — 1-3^0.5 

Fe/Fe0 <0.06 <1.3 0.5lg;3 — — <1.2 — — <0.07 

filling factor . 0.85 0.58 0.74 — — 0.74 — — 0.8 

X^/dof 152/150 147/125 98/85 76/89 129/96 151/121 103/79 113/99 119/116 

^^^Unabsorbed flux in the 0.3-5.0 kcV energy band for the cold (s) and the hot (h) thermal components. The unit is 10^^^ crgscm^^s^^ 



and the site of interaction with the molecular ring (re- 
gion 9); (c) to probe the ISM properties in the northeast- 
ern quadrant, where the surface brightness has its peaks 
(regions 2, 3, 6 and 7), and in the western one, never 
studied before (regions 4, 5 and 8). 

Regions shapes and sizes were determined on the basis of 
statistical and homogeneity criteria so as to have mini- 
mum ^--^5000 counts and small fluctuations of the median 
photon energy (<5%). 

Before performing our analysis, we corrected for the 
vignetting effect assigning a weight to each event (task 
evigweight) and then we extracted source and back- 
ground spectra from the new events list. This proce- 
dure changes the spectral shape of the instrumental back- 
ground, which is not vignetted, so it is necessary to 
extract background spectra in the same regions as the 
observed source data. For fitting we used the sum of 
MOSl and M0S2 spectra, the average of on-axis ARFs, 
weighted for the instruments exposure time, and the 
MOSl canned response matrix ml_110Jm_palLvl.2.rmf. 
The spectra were grouped with a minimum of 25 counts 



per bin and the statistic was used. We first ap- 
plied a single thermal component model, testing both 
for ionization non-equilibrium (VPSHOCK, Borkowski 
et al. 2001) and equihbrium (VMEKAL) cases, and, 
when needed, we added a second thermal component to 
achieve a good fit. We left free to vary O, Ne, Mg and Fe 
abundances of the cold component and Mg, Si, S and Fe 
abundances of the hot component. The two components 
fits with linked abundances give worse results than the 
fits with indipendent abundances. 

Figure shows two representative spectra, extracted 
from regions 1 and 8, with their corresponding best 
fit models and residuals. Best fit parameters of all 
the analyzed regions are presented in Table El with 
the corresponding er rors at the 90% confidence limit 
(jLampton et alJll976(l . 

In the eastern part (regions 1,2,3,6 and 9 in Table EJ 
IC443 thermal emission in the 0.3-5.0 keV band can 
be generally described with a two temperature model: 
the colder (kT~0.3 keV) is charac terized by solar 
abundances l)Anders fc Grevesselll989l) and a ionization 
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timescale t >10^^ cm~'^ s~^, very close to equilibrium 
conditions, the hotter (kT>1.0 keV) by overabundant 
Mg, Si and S and ion equilibrium. On the other side, we 
found the spectra from western regions (number 4,5 and 
8 in Table [SJ to be adequately modeled by a single soft 
component. 

The reduced probability is above 5% for all regions, 
except for region 5 and 9, where we have x^i'=l-3 (1.3%) 
and x^=l-2 (3%). 

4. DISCUSSION 
4.1. Global properties of the X-ray emitting cold plasma 

The soft X-ray emission of IC443 shows high values of 
ionization timescale and a normal chemical composition, 
which suggests an origin from shocked dense circumstel- 
lar or ISM material. 

Previous works about IC443 had left som e ambiguity 
on the plasma ionization state. The work of iPetre et al.l 
l)1988|) . based on Einstein data, was not able to discrimi- 
nate between equilibrium and non equilibrium ionization 
scenarios. Aft erward, the analysis of ROSA T observa- 
tions, made bv lAsaoka &: AschenbachI l)1994j) . described 
the X-ray emission with a two temperatures (~0.3 and 
~1.0 keV) Raymond-Smith model, but the data reso- 
lution did not allow to check for the presence of non- 
equilibrium conditions . On the basis of ASCA data, 
iKawasaki et al.l l)2002t) suggested a high degree of non- 
equilibrium {t^IO^^ cm~'^s) to explain the soft X-ray 
thermal structure of IC443 in its NE quadrant. 
The temperature obtained in the selected regions are 
generally in agre ement with the results of th e "pure 
IC443" regi ons of lAsaoka fc AschenbachI l)1994j) . More- 
over, unlike lKawasakr^t'lLrri)2002() our spectral results 
give a soft X-ray component near the equilibrium condi- 
tion {rr^lO^^ cm~^ s) for all regions except region 7, in 
which we have cm~^ s. The hot component is 

fully equilibrated. 

To infer important values of the gas parameters, such 
as the mean electron density and the total mass of the 
shocked ISM, we assumed a shell-like morphology and 
typical cosmic abundances (ne~1.2 nn) for the soft com- 
ponent. We took into account IC443 double shell struc- 
ture, so we approximated the X-ray emitting volume V 
as two hemispheres with radii Ryi=16' (~7 pc) for the 
eastern side (subshell A in Fig.lSJ and Rb=:26' (^^11 pc) 
for the w estern side (subs hell B in Fig. [SI, following the 
results of lBraun fc StKmj (jll986) and lChevalieil 11999). 

The electron density is thus given by rig ~ EMs/ 
{1.2fsV), where EM^ and fs are the emission measure 
and the soft component filling factor. 
For those regions in which we resolved only a soft 
thermal component, the standard shell geometry was 
adopted: we assumed that the emitting plasma is con- 
fined in a thin shell with a thickness of ~1 pc. To eval- 
uate the filling factors of the other regions we made 
the assumption of pressure equilibrium between the 
soft and the hard component (n,^Tfi'^ nhTh), deriving 
/^i=l-t-(EM^/EM,)(Tft/T,)2 (see e.g. iRocchino et all 
11999). The derived filling factor values are listed in 
Tab. 121 

Fig. 13 shows electron densities and temperatures for 
each spectral region with the corresponding 90% error 
bars. Regions are plotted in order of distance from the 
PWNe. 



Distance from the PWNe (pc) 
1?.5 10.7 9.3 3.6 3.1 5.1 9.5 1-1.7 




6 2 1 9 5 4 

Spectral regions 

Fig. 7. — Electron mean density {top panel) and temperature 
kTs {bottom panel) of the cold X-ray emitting plasma in the nine 
spectral regions indicated on the X-axis. The density values were 
derived from spectral analysis as the square root of the emission 
measure divided by 1.2 times the volume occupied by the cold 
plasma. Subshell A regions (on the left of the vertical dashed line) 
and subshell B regions (on the right of the vertical dashed line) are 
plotted in order of distance from the PWNe. 

Subshell B mean density, with a value of ~1 cm~^, ap- 
pears homogeneous. In subshell A regions nearest to the 
rim (reg 7, 3, 6) there is an indication of a density factor 
two times larger than in the other regions. This is con- 
sistent with an outwardly increasing density profile near 
the rim, which points toward the region of the bright 
very soft shell and optical filaments. 
These results can be explained consistently with a shock 
which is expanding in a homogeneous and low density en- 
vironment in the west, while it has been decelerated by 
the encounter with a density gradient in the east, located 
in a thin region of ~4 pc near the rim. 

We made a rough estimation of the total soft X-ray 
emitting mass taking for each subshell the average of the 
density values, < n >a and < n >b, and assuming that 
the blast wave sheh is 2' thick. We found Mx~30 Mq, 
which is a small fraction of the pr edicted mass of the 
dense molecular cloud (M<10'' M0,'Cornett et al.""197i| 
and the Hi cloud (M'-IO^-IO^ MQ,,Denoyer 1978). This 
is reasonable, also in the light of our results about the 
presence of a density gradient through most of the IC443 
rim (see sect. I4.2(l . since the material emitting in the X- 
rays is just the outskirts of large clouds. 

4.2. Interaction with the environment 

We further investigated the evolutionary stage of the 
X-ray emitting mass in the region between the remnant 
and the molecular clouds, in order to pinpoint possible 
effects of the interaction with these large structures. In 
Fig.|Sl we display contours in the ionization time r vs. 
emission measure EMs space for three regions located in 
sites of putative strong interaction of the remnant with 
the environment, namely region 7 (interaction with the 
neutral clouds), 8 (interaction with ISM) and 9 (interac- 
tion with the southern molecular ring) . We also show the 
post-shock electron density n,, and the age of the emit- 
ting plasma computed using the relation r = n^t, where 
t is the time since the passage of the shock and is 
assumed constant after the shock. The age calculated in 
this way is a very rough estimate and it is more useful 
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for comparing different regions than for the evaluation of 
the true age. 

We stress that for all the regions, the ionization time 
values are compatible with coUisional ionization equilib- 
rium (CIE) at the Scr level. However, we note that region 
7 has a best-fit ionization time which is '--^ 5 times lower 
than 8 and 9 and an age value which is 10 times lower, 
thus opening the possibility that the plasma was more 
recently hit at NE. This is in agreement with a location 
for the SN near the current position of the pulsar wind 
nebula (PWNe), which is also supported by the higher 
values of kTs for the regions near the PWNe (e.g. regions 
1, 5, 9) than for the regions far from it (e.g. regions 6, 
7, 8). According to this interpretation, therefore, the 
regions near the explosion are hotter and more evolved 
than the regions far from it, because they were shocked 
at higher speed long time ago. From this point of view, 
however, it is difficult to explain why the two fronts re- 
gions (7 in NE and 8 in NW) have also different best-fit 
values for ionization time and age, since they are prob- 
ably located at the similar distances from the explosion 
site. 

From the condition pv^=constant, we derived an inde- 
pendent estimate of the evolutionary age of region 7. The 
shock velocity corresponding to plasma temperatures of 
0.3 keV is Vs^^450 km s~^. The density of the X-ray 
emitting plasma that we have calculated is 2.5 cm~^. For 
the Hi cloud we took density and velocity value s from 
the li terature: ri'^lOcm"^ and v~100 km s^^ (Rh o et all 
Unnil)- Optical and X -ray brightness peaks are 23" (2 pc) 
far apart, so we assumed an exponential density distri- 
bution with a scale height of 0.15 pc. The time elapsed 
since the shock front has encountered the cloud is thus 
^2x10"^ yr, in agreement with the hypothesis of a recent 
interaction. 

Alternatively, the low ionization time of region 7 can 
be explained by the fact that the extraction region in- 
cludes plasma at very different densities (e.g. in case of 
strong density gradient). In this case it is possible that 
the VPSHOCK models used in the fit try to reproduce 
the coldest plasma with a low ionization scale. In fact, 
we have seen in Fig. |7| that the density in this region 
increases by a factor of 2.5 in 5 pc (from rig '~1 in reg 1 
to Ue ~2.5 in reg 7), and that there is an additional com- 
pression of a factor of ~4 in 2 pc between the X-ray and 
optical filaments, so this interpretation is realistic. 

A more accurate estimate of the ionization age in both 
regions 7 and 8 would need more data. 

4.3. Absorption from molecular clouds 

The spatial overlap in the image plane between high 
median energy values and molecular emission, shown in 
Fig-EI together with the median/count rate relation sug- 
gest an increase of the foreground absorption in the di- 
rection of CO and H2 emission. 

Using the median photon energy as a spectral indica- 
tor, we selected five representative regions (1, 2, 3, 4 and 
5 in Fig. El to study temperature and column density dis- 
tribution across the FOV. Our best fit results are plotted 
in Fig. ^1 with the corresponding error bars at the 90% 
confidence level. 

They confirm that higher Njj values correspond to higher 
median energies and identify an increasing trend of the 
column density. Moreover, the best fit absorbing col- 
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Fig. 8.— Confidence level ^ contours (68%, 90% and 99%) in 
the EMgtjft-r parameter space for three spectral regions (7, 8 and 
9). On the top X-axis are reported the corresponding post-shock 
electron densities. We show the range of ages which are compatible 
with the values derived from the fit. 

umn of region 1 is significantly higher than the mean 
Galactic absorption toward this portion of sky (N§^' ~ 
6.3x10^^ cm~^, Dickev fc Lockmanl I1990D . indicating 
the presence of additional absorbing material. We es- 
timate that it contributes an equivalent column den- 
sity of ANh'^5x10^^ cm~^. The spatial distribution of 
the obscuring material and the corresponding Njj values 
point to the identific ation with the unperturbed cloud of 
iCornett et all l)1977j) . Our X -ray data strongly suggest 
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Fig. 9. — Two-shells model describing IC443 morphology. The dark ring represents the giant molecular cloud of lCornett et all l|197'i'). 
The blue hemisphere, labelled as subshell A, r epresents the shock front in the eastern region. In the north-east, where it has been confined 
by the encounter with the neutral Hi cloud of 'Denoyei) HQTSf) . the shock front is traced by optical, infrared (2MASS-J energy band) and 
very soft X-ray (0.3-0.5 keV energy band) emission. In the south-eas t the light blue strip indicates the region of interaction with the 
molecular cloud, which instead is bright in the 2MASS-K energy band IRho et alj|200ll) . The orange hemisphere, labelled as subshell B, 
represents the shock front in the western region, where it is expanding in a homogeneous and less dense medium. The ratio of the shells 
radii is 1.6, as derived by our X-ray images. 



that it is located in front of IC443, between the remnant 
and the observer. 

The southern sinuous ridge shows the same E50 val- 
ues of the quiescent cloud, indicating similar hydrogen 
column density. In the median energy map of Fig. |3| 
the southern ridge and the quiescent cloud appear as a 
single complex. We argue that both probably belong to 
the same cloud structure, which has been reached by the 
IC443 shock in several location of the south side (cor - 
responding to clump A to G of 'D ickman et af] 11992V 
while it is mostly unshocked elsewhere, in particular in 
the on-cloud regions ( Fig. El) where there is n o maser or 
H2 infrared emission. iDic kman et alJ l)1992(l suggested 
that the shocked material is in a '^9x7 pc ring tilted 
by ~50° from the line of sight and that it is expanding, 
presumably as a result of the interaction. We expand 
this picture by arguing that the Dickman ring represents 
only the shocked part of a large molecular cloud com- 
plex shaped as a large tilted torus, which is obscuring 
the background SNR X-ray emission. 

In Fig. iniwe show the proposed model for the environ- 
ments of IC443. This model is supported by the com- 
parison (see Fig. between our med ian energy image 
and the contour of lBurton et al.l 1)19881) , which traces the 
boundaries of the shocked gas, and the CO contours in- 
dicating quiescent dense material. While the southern 
and western on-cloud regions are characterized by high 
median energy values (E5o>0.98 keV), which locate the 
molecular ring in the foreground, the northeastern re- 
gions have E5o'^0.9 keV, indicating lower extinction and 



thus, suggesting that the large torus lies in the back- 
ground, as expected in the case of a tilted geometry en- 
compassing the whole SNR. This picture is in agreement 
with kinematic veloci ties measured in clumps A to H by 
IDickman et alJ l|1992|) . and by the m ore recent s tudy of 
clump H and H2 infrared emission of IRho et alJ (|2001). 
In particular, the latter authors detect weak 2MASS-K 
filaments in a region called "NE front", which includes 
clump H and is separated from the bright optical and 
radio NE rim. They argue that the clump H/NE front 
region is part of the molecular cloud interaction, as the 
sinuous southern ridge, and the fact that it is an "off- 
cloud" region according to our Fig. ^] suggests that it 
represents the interaction of the SNR shock with the part 
of the torus lying behind the SNR. 

The comparison between the 0.3-0.5 keV map (Fig.|21l 
and the median energy map (Fig. |2|) shows that the very 
soft shell is not correlated with absorption in the east of 
the remnant, thus indicating that this shell is not entirely 
an effect of absorption. 

The median energy map allowed us to investigate the 
issue of the location of t he nearby SNR G189.6-I-3.3. T his 
remnant was located bv lAsaoka fc AschenbachI l|1994|) in 
front of IC443, causing a dark lane of absorption equiva- 
lent to ANh~2x 10^^ cm~^. This would raise the value of 
the median energy of AE5o^80 eV. Our results does not 
support strong absorption from the shell of G189.6-I-3.3. 
We conclude that either its contribution to Nh should 
be small compared to that of the CO cloud or the SNR 
could be in the background. 
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Fig. 10. — Column density Njj {top panel) and electron tem- 
perature kflT {bottom panel) vs. the median photon energy of the 
indicated spectral regions. The dashed line in the top panel indicate 
the absorption column mean value toward IC443 direction. 

Finally, from the spectral results presented in Table [21 
we argue that the peculiar location of region 5 in the tem- 
perature vs. median energy plot fFig. llOl lower panel) is 
an effect of high temperature rather than column density. 
This is expected, since for Nh< 5x 10^^ cm~^ our median 
energy map is more sensitive to temperature variations 
than for higher Nh values. 

5. CONCLUSIONS 

We presented results from XMM-Newton observations 
of the Galactic SNR IC443. High resolution image in the 
0.3-0.5 keV energy band unveils a partial limb-brightened 
X-ray shell, which resembles on small angular scales the 
optical and radio morphology in the NE quadrant and 
which is d ifferent from the shell structure previously re- 
ported bv lKawasaki et all l)2002D . We concluded that it 



spatially traces the interaction of the shock with a den- 
sity gradient, w hich is pa r t of t he neutral northeastern 
cloud studied bv lDenoveil l)1978|) . 

The spatially resolved spectral analysis of IC443 mea- 
sured large column density variations across the rem- 
nant, and the median photon energy distribution shows 
that they are spatially coincident with molecular emis- 
sion from dense molecular clouds in the FOV. Our re- 
sults confirmed that X-ray e mission is stron g ly abs orbed 
by the molecular cloud of iCornctt ct al. (1977) and 
by the southern sinuous ridge of Burto n et al.. ^988). 
We did not find evidence of absorption from the SNR 
G189.6-h3.3. 

On the basis of our median energy map, we proposed that 
the southern ring and the more extended molecular cloud 
are part of the same structure and that it is tilted re- 
spect to the line of sight, as suggested bv lDickman et alJ 
lll99l . 

We performed spectral analysis in 9 bright regions with 
at least 5000 counts and a 5% maximum variation of the 
median photon energy E50. The X-ray emission from 
our selected regions is well described with CIE or near 
ionization equilibrium models. Our measured X-ray tem- 
peratures and ionization timescales are compatible with 
a density gradient in the NE rim and suggest that the 
site of the SN explosion could be near the PWNe. 
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